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Abstract 

We adapt the methodology proposed by Turfus (2016a) for the analytic pricing of CoCo bonds under 
a Black-Karasinski credit default intensity model to calculate the CVA on a vanilla European equity 
option with a deterministic collateral schedule, incorporating correlation risk between the equity and 
the option writer and a possible jump in the equity price upon default of the option writer. We also 
extend our analysis to cover the case when the collateral is updated to reflect the option value on a 
predetermined schedule. We see that our results are equally applicable to pricing the CVA of a portfolio 
of options. 


1 Introduction 

We look to obtain analytical formulae for the counterparty value adjustment (CVA) for vanilla European put 
and call options on an equity underlying. We take CVA to be defined in the standard manner as described, 
for example, in Hull and White (2012) and Brigo et al. (2013). 

This problem is considered in §8.4 of Brigo et al. (2013). The credit model they used was a structural 
one with assumed correlation between the firm values of the equity issuer and of the option counterparty. 
Zero collateral level was assumed. Results were obtained via Monte Carlo simulation for various levels of 
the correlation, which was seen to have a significant impact on the CVA: a correlation of ±25% can have the 
effect of doubling or halving the CVA. 

Brigo et al. (2013) also consider in §5.3 the CVA on interest rate swap portfolios and swaptions. A 
reduced credit model was assumed for counterparty default risk, based on the Cox, Ingersoll and Ross (1985) 
short rate model. The interest rate model they used was a two-factor Hull and White (1990) short rate 
model. A correlation was assumed between the diffusions of the respective models. Again they calculated 
the CVA via Monte Carlo simulation. They also report in §4.4 analytical results for the swap portfolio case, 
based on swaption modelling and a “drift freezing” approximation, but only under the assumption of zero 
rates-credit correlation. 

We work with a reduced credit model for the option counterparty, supposing that the default intensity 
A t is governed by a Black and Karasinski (1991) short rate model. The equity price is taken to be given 
by a jump-diffusion process, with a downward jump of a fixed relative amount occurring at the time of 
default. The diffusive processes are potentially correlated. We assume from the outset that interest rates 
are deterministic, this being a common assumption in equity option pricing. 

The modelling of the equity and credit default intensity underlyings is presented in §2.1. We consider 
first the result for a call option, presenting in §2.2 the payoff expression and that for the loss on default. 
The equation governing the dynamics of the CVA value is derived in §2.3 supposing in the first instance that 
the collateral level follows a predetermined schedule over the life of the option. The solution for the CVA 
is then obtained using a perturbation expansion approach, as described in §3. It is shown in §4 how the 
approach can be extended straightforwardly to put options and, with a bit more effort, to obtaining similar 
expressions in the case when the collateral value is updated at regular intervals to reflect the current value 
of the option. It is also elucidated how a portfolio of options under a netting agreement can be handled 
analogously. Finally some directions for future work are presented in the closing section §5 
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2 Mathematical Model 


2.1 Basic Framework 

We take the credit default intensity to be governed by a Black-Karasinski short rate model. We shall find it 
convenient to work not with X t directly, but with an auxiliary process y t satisfying the following canonical 
Ornstein-Uhlenbeck process: 

dy t = —ay t dt + a y (t) dW ] (1) 

with yo = 0, where a is a positive constant specifying the mean reversion rate, 1 a y (t) is a bounded positive 
L 2 function and dW^ is a Brownian motion for t £ D m := [0, min{r, T m }] where r is the time of default 
of the option writer and T m is the maturity date of the option under consideration. This equation is well 
known to have a unique strong solution subject to these assumptions. This auxiliary variable is related to 
the credit default intensity X t by 

Xt = (X(t) + X* (t))£(y t ), (2) 

where A (t) is the credit spread defined by 


E 


D - J* \ s ds 


= g- fo A ( s ) ds 


(3) 


and £ (X t ) := exp (X t — 1 [X] t ) is a stochastic exponential with [X] t the quadratic variation of a process X t . 
The required form of the configurable function X*(t) is determined by calibration of the model to satisfy the 
no-arbitrage condition set out below. We see immediately from differentiating Eq. (3) and setting t = 0 that 
Ao = A(0) whence we deduce A*(0) = 0. 

For the equity process we propose: 

dS t j (r(t) — q(t) — kX t ) dt + cri(t) dW? + kdn t if t < r, 

S t . \{f(t)-q(t))dt + a 2 {t) dW? if t > r, 

where St is the equity price, q(t) its expected (continuous) dividend rate, f{t) the instantaneous forward rate 
of interest, dW? a Brownian motion satisfying 

corr(W t 1 , Wt) = pxs (5) 


and nt a Cox process with intensity A t giving rise to an equity price jump of size k with — 1 < k < 0, 
contingent on default. We denote by T w the filtration generated by and Wf and by T J that generated 
by nt . We further define a combined filtration T = T w V T J and formally assume the hypothesis H of 
Elders and Schonbucher (2006) under which T- and -martingales can be considered equivalent. 

We also define a pre-default equity term variance 


rt2 

h(ti,t 2 ) ~ / al(u)du. 


( 6 ) 


If a default occurs at time v with < v < t, 2 , the equity term variance (ignoring the jump impact) will 
reflect the different equity dynamics in Eq. (8) pre- and post-default and be given by 


fV nt2 

I 2 (ti,v,t 2 ) := / al(u)du+ / cr 2 (u)du. 

J t\ J V 


(7) 


When we have t\ = v in Eq. (7), we shall write, for convenience, I 2 (v,t 2 ) := I 2 (v,v,t 2 ). 

Tt is a straightforward matter to loosen this assumption and specify instead that a(t) be a bounded positive L 1 function, 
whereupon the analysis below goes through effectively replacing a(v — u ) throughout by /J a(s)ds. 
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As with the credit default intensity, it is convenient to express the equity price in terms of a new auxiliary 
process Xt, conditional on default at time t = v through 

5 _ "s' 110 ' 0 if t<V, 

1 [ (1 + k)F 2 {v, if t > v, 


where 


Fi(t) := 

Soe-I o( r ( s ) _ 9( s ) _feA ( s )) ds ; 

(9) 

F 2 (v,t) := 

Fi(u}e^ (5;(s) " 9(s))ds , 

(10) 

noting that by construction Xq = 0. We define also for convenience 

(11) 

M 2 (x, v, t) 


(12) 

and risk-free and risky discount functions 

D(ti,t 2 ) 


(13) 

B(ti,t 2 ) 


(14) 


To allow the model to be calibrated we suppose that a term structure of risky bond prices B(0,ti) is known 
for t\ < t -2 < ■ ■ ■ < t n . The required no-arbitrage condition is then that 


D(0,ti)E 


e ~ Jo * ds 


B{0,t z ), i = 1,2, 


(15) 


Eqs. (3) and (15) taken together with an interpolation scheme for times between the nodes fj allow in 
principle both A (t) and A *(t) to be determined for t < t n . 


2.2 Payoff Definitions 

For the CVA calculation, we have two payoffs to consider. In the first instance there is the option payoff 
at at its maturity date, T m say, conditional on no default. Secondly, there is a loss on default at time r if 
r < T m , driven by the value of the option payoff. The latter gives rise to the CVA. We shall consider these 
payoffs in turn. 

In practice we shall need the value of the option not at the initial time 0 but at the default time r = v, 
say. To that end the calculation of the option value will be based on the second line of Eq. (4). Let us suppose 
that x v = £ and y v = rj. The equity value pre-jump at v~ will be determined by the first line of Eq. (8), 
whence that immediately after the jump has occurred will be given by (1 + v) = (1 + k)AI- 2 (f, v , v). 

On this basis the call option price will be given by the standard Black formula, in the present notation: 

fm(Z,v) = D(v,T m ) ((1 + k)M 2 (Z,v,T m )N(d - KN(d 2 (f,v,T m ))) , (16) 


where 


ln((l + k)M 2 {x,t,u)) — In K 

’ ( 

d 2 {x, t , u) + \J I 2 (t, u) (18) 

and iV(-) as usual represents a standard cumulative normal distribution. 

In considering the CVA payoff at time v, suppose that a predetermined schedule governs the collateral 
level and denote this by C(t). The CVA loss calculated at v will be given by 1 — R times the difference 


d 2 (x,t,u) := 

di(x,t,u) := 
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between the option value and this collateral level if positive, zero otherwise, where R is the assumed recovery 
level of the counterparty. Thus the payoff function for default at time r = v is given by 


P m (£,v) = (1 - fi)max{/ ra ((,w) - C'(v), 0} . 

The counterparty value adjustment at time t > 0 is then defined in our notation to be 


(19) 


CVA(i) = (1 - R) / D(t,v)E[l r > v P m {x v ,v)\F t ]dP{v) 


= (1 -R) / m D(t,v)E\ \ v e-K x ’ ds P m {x v ,v) 


•F, 


w 


dv, 


( 20 ) 


where P(-) is the probability distribution of the time of default events driven by the Cox process n t defined 
above and in deriving the last line we have relied upon the result of Elders and Schonbucher (2006). 


2.3 Main Equation 

We next suppose that the process CVA(t) can be expressed in terms of the new (stochastic) variables Xt and 
yt as h(xt,yt,t). We are interested in calculating CVA(O) = h(0, 0, 0). Applying the well-known Feynman- 
Kac method to Eq. (20) subject to Eqs. (1), (4) and (5), we infer that the function h(x,y,t) representing 
the CVA is governed by the following backward diffusion equation: 


C[h(x,y,t)\ 


(A (i) + AA (y, t)) P m {x, t) + AA (y, t) 



( 21 ) 


for t £ D rn with final condition that h(x,y,T m ) = 0, where £[•] is a forced diffusion operator given by 

£H : = I - ay Wy + l |?) - (F(1) + J{t)) (22) 


AA (y,t) := (A(t) + \’(t))£ y ( m ) - A ((), (23) 

where £ y (y t ) for y £ K is a shorthand for £(yt)\y t =y Our task is then to solve Eq. (21) to obtain 
CVA(0) = h(0, 0, 0). 


3 Perturbation Analysis 


In the absence of exact closed form solutions to Eq. (21), we seek an approximate solution for CVA (t) under 
the low default intensity assumption proposed by Turfus (2016b), defining an asymptotically small parameter 
e by 

1 


aT, 


X(t)dt. 


m J 0 


We then define a scaled forward intensity A (t) by 

A (t) := e _1 A (t) 

and further define 


taking these to be 0(1) as 


AA (y,t) = e AA (y,t) 

0. Turfus (2016b) shows that we can write 

X*(t) = e 2 X* 2 (t) + 0(e 3 ). 


(24) 

(25) 
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The form of A |(1) is not required for the present purposes so need not concern us further here. 

To solve Eq. (21) in it asymptotic representation we pose a perturbation expansion 

h(x,y,t ) = ehi(x,y,t) + e 2 h 2 (x,y,t) + 0(e 3 ). (26) 

Substituting the above expansion into Eq. (21), solving up to 0(e), setting 1 = 0 and reverting to unsealed 
notation, we obtain the following first order accurate result: 

Theorem 3.1 The CVA on a European call option on an equity underlying can be estimated under our 
modelling assumptions as follows: 

CVA(0) = (1 -R) [ " B(0,v)D(v,T m )X(v)(l + k)F 2 (v,T m )e I ^°’ v) ^+(v,T m ) dv 
Jo 

-(1-1?) [ B(0,v)X(v) (D(v,T m )K^(v,T m ) + C(v)N(a 2 (I p (0,v),0,v))) dv 

Jo 

+ 0(e 2 ) (27) 


I p (ti,t 2 ) := P xs J t e Cy(u)cTi(u)dv 


1 pf(v,w) := N 2 (±a,i(Ip( 0,v),0,v),±bi(I p (0,v),0,v,w),R(0,v,w)), 1 = 1,2, 


a 2 (x,t,v) := 


x-C(v) 


R(t,v,w ) : = 


v < w, 


ai(x,t,v) := a 2 (x,t,v) + \/li(t,v), (31) 

ln((l + k)M 2 (x, v, w)) — In K 

b 2 (x,t,v,w) := , (32) 

V h(t,v,w ) 

bi(x,t,v,w) := b 2 (x,t,v,w ) + \J I 2 (t, v , w) , (33) 

R(t,v,w) := J Jl ^’ V - , v<w, (34) 

V I 2 {t,v,w) 

C(v) ■■= sup{^|P m (C,'c) = 0}, (35) 

and N 2 (xi, x 2 ; p) is a standard bivariate cumulative Gaussian distribution with correlation p. 

Proof. For the proof of Theorem 3.1, see Appendix A. Note that in terms of our previous notation 

bi(x, v , v, w) = di( x, v, w), 1 = 1,2. □ 

The corresponding result for the CVA on a put option is obtained similarly. 

Theorem 3.2 The CVA on a European put option on an equity underlying can be estimated under our 
modelling assumptions as follows: 

CVA( 0) = (1 - R) [ m B{0, v)X(v) (D(v, T m )K (v, T m ) - C(v)N(-a 2 (I p ( 0, t>), 0, v))) dv 
Jo 

-(1-1?) f ”s(0,^A(t;)(l + fc)lA(u,T m )E 2 ( ? ;,T m )e 7 ^ 0 ’ t ’VrKr m )du + O(e 2 ). (36) 

-'O 

The notation is as defined for Theorem 3.1, except in that the option price must now be reinterpreted as 

fm&v) = D(v,T m ) (. KN(-d 2 (Z, V ,T m )) - (1 + fc)M 2 (^,t;,T m )lV(-d 1 (^r;,T ro ))) (37) 

in the definition of P m (£,v), rather than as given for the call option in Eq. (16). 2 

2 We must also reinterpret £*(u) as inf rather than sup contingent on this reinterpretation of notation. 
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Proof. The proof follows closely that of Theorem 3.1 in Appendix A, making use of the formula for a 
compound call option on a put option provided in Geske (1979). □ 

In interpreting the significance of these results, we note that the main impact of correlation will be 
through the appearance of the exponential involving / p (-) multiplying the equity forward. If the equity and 
credit intensity are assumed negatively correlated so that l p (0,v) < 0, the equity price at default will be less 
than otherwise assumed. The value of the call option at default will in this way be reduced, whence the risk 
is inferred to be “right-way” (and correspondingly “wrong- way” for the put). Concerning the impact of the 
(downwards) jump at default, the main consequence of this will be further to decrease the value of the equity 
at default so will impact in the same direction as the assumed negative correlation. For short-term options, 
the jump effect will in likelihood be the more important: the impact on world stock markets from the failure 
of Lehman Brothers was around 3-4% downwards. But as the impact of correlated diffusion on the option 
price will accumulate over time with, by our assumption, no compensating adjustment made to the collateral 
level, for longer maturity options, this effect will ultimately exceed the jump impact and become the more 
important. 


4 Extensions of Main Result 

4.1 Dynamic Collateral 

It is of course often the case in practice that collateral requirements are not fixed at the onset of the trade 
but change to reflect the level of exposure. In this section we look to see how the results of Theorems 3.1 
and 3.2 can be extended to take account of a policy whereby at scheduled dates tj, j = 0, 1, . . . , n — 1 the 
collateral is updated to reflect the current exposure level associated with the option, viz. its PV. Clearly in 
this situation, each exposure period [tj,tj+ 1 ) can be considered separately, with the compound option taken 
to be forward-starting at tj with strike given by the value at time tj of the option, conditional on survival 
of the option writer. 3 To that end we need to take the market view of the option value at time t > tj 
conditional on survival to time tj. We argue that this can be achieved by considering the equity dynamics 
to be given by the first line of Eq. (4) until tj and by the second line thereafter. Thus we propose 

S t = F 2 (t j ,t)e Xt ^ l2{0 ' t ^ t \ t>tj> 0 

By straightforward calculation, we infer the market value of the option as of time tj conditional on x t . = x 
and t > tj is 

Cj(x) = Ditj^^&tj^Nidl&tj,^)) - KN(d* 2 (x,tj,T m ))), (38) 

where 


d 2 (x,t,u) := 

In M 2 {x, t, u ) — In K 

(39) 

\J h{t,u ) 

d\{x,t,u) := 

d 2 (x,t,u) + I 2 {t, u) 

(40) 

The payoff function we must now consider conditional on Xt : Xj for j = 0, 1, . . . , n — 1 becomes 


:= (1 - R) max {/ m (£, v) 

+ 

W"i 

1 

tcf 

1 

1 

S 

o' 

B 

(41) 


with H(-) the Heaviside step function. In this way we obtain: 


3 We have ignored the possibility that to ^ o, viz. the situation where the collateral for the first period is already known 
rather than being determined dynamically. However this is the case already considered in Theorem 3.1 so can be dealt with in 
the manner described there, with the caveat that ti rather than T m must be used as the upper limit in the integrals. 
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Theorem 4.1 The CVA on a European call option on an equity underlying associated with default in 
a forward time period [tj,tj+ 1 ) where the collateral is reset at tj can be estimated under our modelling 
assumptions and the approximation Eq. (55) as follows: 

CVA & (0) « (1 -R) f 3+1 B( 0, v)D(v, T m ) X(v)(l + k)F 2 (v, T m )e / ^°^# + (^(0, v), v, T m ) dv 

J tj 

- (1 - R)K [ 3+1 B(0,v)D(v,T m )X(v)^ j)+ (I p (0,v),v,T m )dv 

Jtj 

-(1-1?) f X(v)B(0,v)D(v,T m )N (cL$\l p {tj,v),t j ,v) s j dv 

Jtj 

(e /p(0 ^ ) F 2 (t j ,T m )N(b[ j \l p (0,tj),0,T m )) - KN(b { 2 j) (I p (0,tj), 0, T m ))) , (42) 

while the corresponding result for the CVA on a European put option is 

CVA^X 0) « (1 - R)K [ 3+1 B{0,v)D(v,T m )X{v)^ j) -{I p (0,v),v,T m )dv 

J tj 

- (1 - R) f J+1 B(0,v)D(v,T m )X(v)e I ^ v Xl + k)F 2 (v,T m )tp[ j) -(I p (0,v),v,T m )dv 

Jtj 

-(1-1?) f X(v)B(0, v)D(v, T m )N (-a^Xl p (tj,v),tj,v)^ dv 

Jtj 

(^KN(—bXXl p (0,tj), 0, T m )) — e Ip( -°’ tj ^F2(tj,T m )N(—biXl p (0, tj), 0,T m ))j , (43) 

where 



■= N 2 {±a[ J Xl p (0,v),0,v),±b ( i J Xl p (0,v),0,v,w), R(0,v,w)), i = 1,2, 

(44) 

a { 2 J \x,t,v) 

x-€j(v) 

' n? — t. r i 

(45) 


Vh{t,v) 


a^\x, t, v) 

:= a 2 (x,t,v) + y/h{t,v), 

(46) 

b { 2 j) (x,t,w) 

In M 2 (x, tj , w) — In K 
' / T / \ 5 

(47) 


V h{t,tj,w ) 


bi\x,t,w) 

:= b 2 (x,t,w) + Jl 2 {t,tj,w) 

(48) 


and the definition of £* (y) is given in Appendix B below. 

Proof. For the proof of Theorem 4.1, see Appendix B. □ 

The total CVA for the option is then simply the sum over j of all contributions of the form of Eq. (42) or 
(43). Cancellation between terms means that the impact of correlation is considerably reduced in comparison 
with the deterministic collateral case considered in Theorems 3.1 and 3.2 above. Basically the impact of 
correlation on CVA only has the chance to build up between collateral readjustment dates, when it is reset 
to zero. However, the impact of the jump is expected to stay as strong as previously, so will tend always be 
the greater of the two effects. 

4.2 Portfolio of Options 

The extension of the above result to a portfolio of options on a single underlying, potentially combining calls 
or puts with various strikes and maturities is not difficult. We assume here that the collateral for all options 
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is updated on a single schedule. We further assume that a collateral adjustment occurs on the expiry dates 
of all the options. Rather than using the call option payoff of max{/ TO (£, v) — f m (xj,tj), 0} in Eq. (41), we 
must use a sum over all the options in the portfolio, namely 

max | (fm (£’ v ) ~ C f Oj )) i 0 1 

where the superscript (*) denotes the ith option in the portfolio. Reinterpreting £*(xj,v) = Xj + £*(v) as 
relating to this aggregated payoff function for a given v £ [tj, ij+i), the analysis goes through identically to 
obtain the CVA of each call option separately using the same £* (v) in all cases; similarly for a portfolio of 
puts. If there are both calls and puts, it is expected to be the case that the payoff will be in the money if £ 
is high enough or low enough. So there will be two critical values, each giving a contribution to CVA^(O) 
when its value is exceeded. There is nothing to prevent the inclusion either of both long and short positions 
in the portfolio, provided it is the case that collateral is posted to the counterparty (as usually would happen) 
should the total value of the portfolio ever become negative. 

5 Future Work 

It, has not yet been possible to assess quantitatively the relative importance of the correlation and jump 
correction terms proposed above for the various cases considered, nor to demonstrate the accuracy of the 
assumed first order expansions or the impact of the approximation Eq. (55). This is the subject of ongoing 
research work by the author. 

It should also be possible to extend these results fairly straightforwardly to other European equity options 
such as digitals. Similarly the equity underlying could be swapped for an FX rate, an inflation rate or a 
commodity spot price. 


A Proof of Theorem 3.1 


We derive Eq. (27), solving Eq. (21) in the standard manner by successive levels of approximation. At first 
order we must solve 


£[hi (x,y,t)\ = -X(t)£ y (y t )P m (x,t) 


for t £ D m with final condition that h\(x,T m ) = 0. This can be achieved by means of the following readily 
obtainable Green’s function for the diffusion operator £[■]: 


d 2 

G(x,y,t;£,v,v) = B(t,v)-Q^N 2 


ye “ 1 ” — 77 


x-j 

\/h{t,v)' y/Iy(t,v) 


; \p(t,v ) 


0 < t < v, 


(49) 


where B(t,v) is the risky discount factor defined in Eq. (14) above, 


P(t,v) 
1 y (G > 1 2 ) 


I P (t,v) 


^/h{t,v)ly(t,vy 

pt2 


e -2a(t,-u) £r 2( u ) du 


(50) 

(51) 
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and Ii(ti,t 2 ) and I p {t\,t 2 ) are as defined in Eqs. (6) and (28) above, respectively. Making use of Eq. (49) 
we obtain 


hi{x,y,t) 


(1 - R) 




- (7(r;)) f £r,(y v )G(x,y,t;£,ri,v)dridtdv 


(1 - R)£ y (y t ) J \{v)B(t,v) 

( e x+/ p (t,«)- i/, (0,t) (1 + k)D ( v ^ Tm )F 2 (v, T m ) 

N 2 {a l (x + I p (t, v),t, v), b\(x + I p {t, v),t, v, T m ), R(t, v, T m )) 

- D(v, T m )KN 2 {a 2 {x + I p {t , v), t, v),b 2 (x + I p {t, v), t , v, T m ), R(t, v, T m )) 


— C(v)N(a 2 (x + I p (t,v),t,v)) S jdv. 


(52) 


Here in performing the integration w.r.t. £ we have followed the analogous calculation by Geske (1979) of 
the price of a compound call option on a call option. 

This process can be continued to second order to obtain h 2 (x,y,t) but we do not consider such higher 
terms here. Setting x = y = t = 0 in Eq. (52) and reverting to unsealed notation, we obtain Eq. (27). This 
concludes the proof of the theorem. □ 


B Proof of Theorem 4.1 


We seek to derive Eq. (42) as the solution to Eq. (21) with the exposure function given by Eq. (41) for 
v € [tj,tj+ 1 ). Let us take h^\xt,yt,t) for t < tj as the value at time t of the CVA contribution associated 
with defaults in the interval [tj,tj+i) and suppose it can be written 

h {j) {x, y, t) = eh[ 3 \x, y,t) + e 2 h ( 2 J \x, y, t) + 0(e 3 ) (53) 

for 0(1) functions The final condition satisfied is h^\x,y,tj + 1 ) = 0. It is convenient to calculate the 

desired (first order) result (0,0,0) in two stages. First, estimate h[ 3 \x,y,tj): as the underlying is not 
then a forward-starting option, the calculation is almost exactly as was performed before. Next we use the 
tower property of conditional expectations to assert that the CVA at time t < tj associated with the future 
period [tj, tj+i) is given by considering the same Eq. (21), but with the final condition that h[ 3 \x,y,tj) is 
as previously determined. The desired result is then obtained by setting t = 0. 

Following the approach set out in Appendix A, we find from the first order calculation that 

hi\x,y,tj) = (1 - R) [ A (v) f (fm{£,v)-C j (x)\ ( £ ri (y v )G(x,y,t j ;^,r],v)dr]d^dv, (54) 

J tj j£*(x,v) ^ ' J— oo 

where 

Cj{x,v) := inf{£ | f m {£,i v) - Cj{x) > 0}. 

By its construction, we expect £*{x,v) to be approximately of the form 

£j(x, v )=x + £*(0,v). (55) 

We will henceforth make this approximation, which simplifies the calculation considerably and adopt the 
convenient notation that £*( v) := ^*(0,f). On this basis we obtain 


h[ j) (x,y,tj) 


(1 - R)£ v (y tj 


hl A (v)B(tj,v) f°° 

\/l\ (tj , l)) J£*(x,v) 


N' 


( x + I P {tj,v)-( \ 

V J 


fm(£,,v)d£dv 


(1 - R)Cj{x)£y{y tj ) 



\{v)B(tj,v)N 


V \Aife.f) J 


dv 


(56) 

(57) 
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Applying our Green’s function straightforwardly to a payoff of h^\x,y,tj) at final time tj we obtain 


/>oo />oo 

hi\x,y,t)= / G(x,y,t;£,r),t j )h ( ? ) (t,r},t j )dr)d4 

J — oo J — oo 

= (1 - R)Sy(yt) j t J+1 A (v)B(t, v)D(v, T m ) (>+C(M)-|/i(o,t) (1 + ^ 

N 2 (a[ J \l p (t, v),t, v ), 6i(a; + I p (t, v),t, v, T m ), R(t, v, T m )) 

- KN 2 (a%\l p (t, v),t, v), b 2 (x + I p (t, v), t, v, T m ), R{t, v, T m )))dv 


(1 - R)£ v (y t )J 3+1 A (v)B{t,v)D(v,T m )N dv ^+G(*,b)-|A(o ,t) F ^ tjyTm ) 


N(b[ j) (x + I p (t, tj),t, T m )) - KN(b^’(x + I p (t, tj), t, T m )) 


U)i 


(58) 


for t < tj. Setting x = y = t = 0 and reverting to unsealed notation, we obtain Eq. (42). This concludes the 
proof of the theorem. □ 
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